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FOREWORD 


This  report  documents  an  analytical  effort  for  the  estimate  of  the  effect 
of  wall  heating  on  the  transition  location  of  the  boundary  layer  in  a tube. 
This  effort  complements  the  DARPA-sponsored  flow-tube  experiment  being 
conducted  at  the  Colorado  State  University.  The  findings  reported  herein 
clearly  identify  several  critical  aspects  of  the  interpretation  of  the 
experimental  results  acquired  in  the  flow-tube  facility.  Several  rather 
favorable  comparisons,  in  contrast  to  the  apparent  discrepancies  cited  in 
the  past,  with  the  test  data  are  obtained  when  the  present  findings  are 
taken  into  consideration. 

This  work  was  sponsored  by  DARPA/TTO  under  ONR  Contract  N00014-77-C-0005 
( P0002 ) . 
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1.  INTRODUCTION 


In  recent  years,  several  experimental  efforts  have  focused  on  the  use 
of  surface  heating  as  a means  of  stabilizing  the  laminar  boundary  layer 
adjoining  a body  moving  in  water.  In  order  to  obtain  data  in  a more 
controlled,  oteady-state  environment,  Autonetics/Rockwell  International 
has  initiated  a series  of  tests  at  the  Colorado  State  University  flow- 
tube  facility.  The  primary  results  obtained  from  this  series  of  ex- 
periments are  summarized  in  Figure  1 (Barker,  1977).  A substantial 
increase  in  the  transition  Reynolds  number  is  evident  when  the  surface 
overheat,  AT,  is  increased  from  0 to  5°C.  Above  5°C,  the  effect 
of  heating  appears  to  be  less  consequential  and  the  transition  Reynolds 
number  approaches  an  apparent  upper  limit  of  R^  = 42x10^.  The  theor- 
etical predictions  of  Wazzan  et  al  (1968,  1970)  suggested  that  transi- 
tion Reynolds  numbers  of  up  to  200x16^  may  be  possible  in  a zero 
pressure  gradient  boundary  layer  with  surface  heating.  Dashed  lines 
in  Figure  1 show  the  results  from  those  predictions  for  a zero  pressure 
gradient  boundary  layer  (6=0)  and  for  a boundary  layer  subjected  to  a 
favorable  pressure  gradient  (8=0.07).  The  rather  significant  disparity 
between  the  predictions  and  the  experimental  results  has  been  the 
primary  cause  for  numerous  experimental  modifications. 

In  a recent  study,  Tzou  et  al  (1977)  indicated  that  the  noted  discrepancy 
between  the  theoretical  predictions  and  the  experimental  results  might  be 
attributable  to  the  geometric  differences  between  a flow  tube  and  a flat 
plate.  These  conjectures,  however,  must  be  further  substantiated  by 
comparison  of  the  experimental  findings  with  more  accurate  analytical 
predictions  of  the  transition  location  in  a flow  tube.  The  main  features 
unique  to  a flow-tube  system,  such  as  the  effects  of  transverse  surface 
curvature,  upstream  contraction  section  and  boundary-layer  displacement 
thickness,  should  be  included  in  these  "exact"  computations.  These  main 
features  unique  to  a flow-tube  are  discussed  in  Section  2.  Results  of 
the  calculations  for  a heated  tube  with  an  unheated  extension,  using  this 
improved  theoretical  prediction  method,  are  presented  in  Section  3.  A 
summary  and  several  conclusions  based  on  the  present  results  are  in 
Section  4. 


Overheat 


Wall  Overheat  AT(°C) 

Figure  1.  Boundary  Layer  Transition  Data  and  Predictions  for  the  Colorado 
Flow  Tube  Experiment  (Barker,  1976,  1977) 
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2.  MAIN  FEATURES  UNIQUE  TO  A FLOW  TUBE 

A complete  simulation  of  a hydrodynamic  phenomenon  requires  that  the  two 
systems  be  geometrically,  kinematically,  and  dynamically  similar.  Because 
of  physical  constraints,  it  is  impossible  to  achieve  an  exact  simulation 
of  the  flow  field  around  a given  external  axisymmetric  body  with  an  internal 
flow- tube.  In  order  to  be  able  to  utilize  the  experimental  results  obtained 
from  a flow  tube,  a thorough  understanding  of  the  flow  characteristics 
within  a tube  is  essential.  Several  key  elements  unique  to  a flow-tube 
system  have  been  identified  (Tzou  et.  al . , 1977).  These  elements  need  to 
be  included  for  more  accurate  computations  and  they  are  briefly  discussed 
below. 

2.1  Transverse  Surface  Curvature 

Generally,  the  effect  of  transverse  surface  curvature  is  important  only  if 
the  boundary  layer  thickness,  <5,  is  not  negligibly  small  as  compared  to  the  body 
radius  r . Furthermore,  the  effect  is  different  for  internal  and  external 
flows  as  is  discussed  below. 

In  the  following  sketch  "r"  represents  the  distance  measured  from  the  body 

axis  and  is  related  to  "r0",  the  local  radius  of  the  body  surface  (i.e., 

y = 0)  through  r = rQ  + y cosa  for  an  external  flow  and  r = rQ  - y cosa  for 

an  internal  flow,  where  a is  the  angle  that  the  surface  makes  with  the  body 

axis.  A transverse  curvature  term,  t = ^ cos  a , is  introduced.  Then,  for 

r 

an  external  flow: 

= 1 + t , and  (1) 

o 

for  an  internal  flow: 


(2) 


EXTERNAL  FLOW  INTERNAL  FLOW 


In  order  to  investigate  the  quantitative  effect  of  this  t term  in  the 
boundary- layer  equations  for  an  internal  flow,  the  TAPS  Code  (Gentry  & Wazzan, 
1976)  is  modified  to  include  the  effect.  The  effect  of  transverse  curvature 
on  the  evolution  of  a boundary-layer  is  illustrated  by  considering  the  Rockwell/ 
Autonetics  flow-tube  which  is  circular  in  cross  section  with  a diameter  of 
10.16  cm  and  is  6.10  m long.  The  flow  conditions  used  tor  *his  demonstration  are: 


freestream  speed  UQ  = 6.10m/sec 

freestream  temperature  T^  = 10°C 

wall  overheat  AT  = Tw-Tro  ='5.56°C. 

In  order  to  isolate  the  effects  of  the  transverse  surface  curvature,  the 
pressure  coefficient,  Cp(x),  is  assumed  to  be  zero  for  both  the  internal 
flow  as  well  as  the  equivalent  external  flow.  The  computed  boundary-layer 
characteristics  are  presented  in  the  following  table. 
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Note  that  the  effect  of  flow  geometry  on  the  boundary-layer  thickness  is 
quite  small  even  when  the  ratio  of  the  boundary-layer  thickness,  6,  to 
the  radius  of  the  tube,  rQ,  is  about  0.11  at  this  location.  The  effect  on 
the  momentum  thickness,  6,  and  the  shape  factor,  H,  is,  however,  significant. 
It  is  also  interesting  to  reveal  that  the  transverse  surface  curvature  causes 
an  internal  flow  to  have  a higher  value  of  shape  factor  (H  = 2.72)  than  its 
external  flow  counterpart  (H  = 2.40). 

2.2  Upstream  Contraction  Section 

In  the  cited  flow-tube  experiments,  a contraction  section  exists  upstream 
of  the  flow  tube  in  order  to  reduce  the  60.96cm  diameter  supply  pipe  to 
the  10.16cm  diameter  test  section.  As  shown  in  Figure  2,  the  boundary-layer 
displacement  thickness  in  the  contraction  section  varies  with  the  flow 
conditions.  Therefore,  for  an  accurate  calculation  of  the  flow  characteris- 
tics in  the  tube,  calculations  should  begin  upstream  of  the  contraction 
section. 

2.3  Boundary  Layer  Displacement  Thickness 

Because  of  the  physical  constraint  of  the  boundaries  in  an  internal  flow,  a 
much  stronger  effect  on  the  freestream  flow  is  expected  than  for  its  external 
flow  counterpart.  Since  the  flow  rate  Q is  constant  along  the  flow  tube  but 
the  effective  cross-sectional  area,  itr2  = Tr(ro-6*)z,  decreases  in  the  down- 


stream direction  (see  sketch  below)  as  a result  of  boundary-layer  growth,  the 
"freestream"  speed  must  change  according  to: 


We  note  that,  for  the  present  purposes,  simple  one-dimensional  flow  relations 
are  used.  The  effect  of  this  speed  change  on  the  pressure  coefficient,  Cp, 
can  be  expressed  as 


Figure  3 shows  the  variation  of  Cp  along  the  Autonetics  flow-tube  for  two 

different  conditions.  Without  applying  the  displacement  thickness  correction 

on  C , it  is  customary  to  assume  C = 0 throughout  the  flow-tube.  With  such 
P P 

a correction,  Cp  at  the  end  of  a 6.10m  flow  tube  is  -0.13  when  U0  = 6.71m/sec, 
AT  = 16.67°C  and  Cp  is  -0.27  when  UQ  = 1.83m/sec,  AT  = 0°C.  It  is  significant 
to  note  that  the  general  assumption  of  either  a zero  pressure  gradient  or 
a constant  pressure  gradient  flow  within  the  tube  is  definitely  incorrect. 

Since  the  displacement  thickness,  6*,  is  a consequence  of  the  particular 
pressure  distribution  in  the  flow-tube,  an  accurate  estimation  of  6*  can  be 
obtained  only  by  an  iterative  procedure.  To  initiate  such  an  iterative 
calculation,  it  is  assumed  that  6*  can  be  first  estimated  by  (See  Figure  4) 


where  XQ  is  the  length  of  the  contraction  section,  A is  the  distance  from 
a virtual  origin  to  the  beginning  of  the  flow  tube  test  section.  It  should 
be  noted  that  this  simplification  was  only  used  to  start  the  computation  and 
that  the  actual  calculation  of  6*(X)  was  used  for  subsequent  iterations. 


Figure  3.  Pressure  Coefficient  Distribution  in  the  Flow  Tube 
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A flow  diagram  showing  the  iterative  procedure  to  account  for  6*  correction 
on  C in  the  boundary  layer  flow  calculation  and  the  subsequent  stability 
analysis  is  described  below  in  conjunction  with  Figure  5.  Note  that  all  three 
features  discussed  in  the  section  are  included  in  these  computations. 

1.  The  computation  requires  input  data  such  as  the  tube  geometry 
X,  rQ  and  the  flow  conditions  U0,  AT,  T«>  and  v. 


2. 


Use  equation  (5)  to  calculate  the  initial  6*  value  by  assuming 
C i =1 . 5 , X0=1.37m  and  A=0.27m 


%) 


= C, 


v(X-X0+A) 


Uo 


3. 


Use  equation  (4)  to  obtain  the  corresponding  variation  of  Cp 
along  the  flow  tube 


CP(D  = 1 ' 


ro  " ^o) 


4,  Perform  Boundary  layer  flow  calculation  to  obtain  calculated 

5.  Compute  new  cp(j+i)  from  the  calculated  6^) 

r r„  1 4 

Cp(i+D  = 1 " 


L - «ti) 


6.  Calculate  the  convergence  parameter 


F C 


£i  " 


p(i+l)  - p(i ) 


'p( i+1 ) 


J 
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7.  (a)  If  the  iteration  i < 3 or  the  convergence  parameter  c.  > 0.005, 

the  procedure  return  to  step  (4)  for  the  next  iteration  (i+1). 

(b)  If  i > 3 and  e-  < 0.005  the  calculation  proceed  to  the  next 
step. 

8.  Perform  stability  analysis  based  on  the  results  of  boundary  layer 
flow  calculation  from  last  iteration. 

9.  Determine  the  Loge  (A/AQ)  for  the  corresponding  flow  conditions 
specified  in  step  (1)  from  the  spatial  amplification  factor  map. 


i 
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RESULTS  - A HEATED  FLOW  TUBE  WITH  AN  UNHEATED  EXTENSION 


All  the  effects  discussed  in  Section  2 have  been  incorporated 
into  the  boundary-layer  component  of  the  TAPS  Code  to  compare  with 
the  experimental  results  obtained  with  one  unheated  extension 
tube  and  a laminar  flow  nozzle  that  are  shown  in  Figure  6.  Four  of 
these  "best"  experimental  flow  conditions  (solid  circles  in 
Figure  6)  are  selected  for  further  analytical  study. 


Flow  Tube  Inlet 
Heated  Tube 
Unheated  Tube 
Ambient  Temperature 


AT(°C) 

Rt  (F 

0.00 

12.0 

2.78 

25.0 

5.56 

40.0 

8.33 

40.3 

1.37m 

6.10m 

1.22m 

Too  = 10°C 


gure  6)  UQ  (m/sec) 


x 

106 

2.06 

X 

106 

4.29 

X 

106 

6.86 

X 

106 

6.89 

Following  the  procedure  described  in  the  flow  diagram  (Figure  5), 
three  iterations  were  required  to  satisfy  the  convergence 
criterion  of  e..  <_  0.005  for  the  four  selective  cases.  As  shown 
in  Table  1 and  also  in  Figures  7 and  8,  the  iteration  of  the 
calculation  converge  very  rapidly.  Therefore,  the  initial  selec- 
tion of  the  coefficients  Ci,  A and  X0  in  equation  5 is  not  so 
critical  in  the  computation.  In  this  study,  same  coefficients 
Ci  = 1.50,  A = 0.27m  and  X0  = 1.37m  are  used  for  initiating 
the  calculations. 

Variations  of  boundary-layer  characteristics  in  the  flow  tube 
are  illustrated  in  Figures  9 through  11. 


OVERHEAT 
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Figure  6.  Experimental  Transition  Reynolds  number  data  showing  scatter 
due  to  water  quality  variation.  (Barker,  1978.) 
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Figure  9 shows  the  variation  of  the  pressure  gradient  parameter  6 
(defined  as  g = , where  £ is  the  transformed  distance  parameter, 

Ug  Qs 

see  Gentry  & Wazzan,  1976)  in  the  flow  tube  for  various  operating  con- 
ditions. We  note  that  the  computed  values  of  g depend  strongly  on  the 
flow  speed  and  increase  considerably  in  the  downstream  direction.  The 
flow  in  the  tube  is  neither  a zero  pressure  gradient  flow  (6  = 0)  nor 
a constant  pressure  gradient  parameter  flow.  Hence  analytical  results 
based  on  either  of  those  assumptions  would  not  be  strictly  legitimate. 
It  is  also  noted  that  parameter  g has  a sudden  jump  at  the  joint  of 
heated  and  unheated  sections. 

Figure  10  shows  the  variations  of  shape  factor  H (H  = 6*/9 ) in  the  flow 
tube.  In  the  cases  of  heated  flow  conditions  (AT  = 2.78°C  and 
AT  = 8.33°C),  the  shape  factors  increase  considerably  in  the  unheated 
extension  tube,  indicating  the  boundary-layer  has  become  more  unstable. 

Figure  11  shows  the  variations  of  the  displacement  thickness  in  the 
flow  tube  for  various  flow  conditions.  As  can  be  seen  in  the  figure, 
a lower  flow  velocity  has  a higher  displacement  thickness  and,  in  turn, 
has  a much  larger  pressure  gradient  parameter  8. 

The  boundary-layer  stability  analyses  have  been  performed  for  the 
selected  flow  conditions.  The  computed  spatial  amplification  factors 
for  various  disturbance  frequencies  are  shown  in  Figures  12  to  15.  A 
summary  of  the  results  is  presented  in  Table  2. 


At  the  experimental  transition  location,  end  of  the  flow  tube  with  one 
extension,  the  corresponding  amplification  factors  vary  from  e11'6  to 
e15'®  with  a mean  value  of  e14'0.  All  of  these  amplification  factors 

9 

are  much  higher  than  the  traditional  "e  " criterion  which  was  proposed 
by  Smith  (Jaffe,  1969). 


12.9 


Figure  12.  Spatial  Amplification  Factors  in  Boundary  Layer  for 
Ub=2.06  m/sec,  AT=0°C,  =10°C 


U0=6.86  m/sec,  AT=5.56°C 
T =10°C 


Figure  14.  Spatial  Amplification  Factors  in  Boundary  Layer  for 
U0r6.86  m/sec,  AT=5.56°C,  T =10°C 
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TABLE  2 

Amplification  Factors  At  The  End  Of  One  Extension  Tube 


Non-Dimensional 
Frequency 
u x 105 

AT=0°C 

Uo=2.06m/sec 

AT=2. 78°C 

U0 =4. 29m/ sec 

AT=5.56°C 

U0=6.86m/sec 

AT=8.33°C 

U0=6.89m/sec 

2.25 

7.63 

2.00 

8.50 

1.75 

9.41 

1.50 

10.83 

6.43 

1.25 

* 12.91 

7.70 

3.11 

0.07 

1.00 

11.08 

9.50 

4.37 

0.49 

0.90 

9.14 

10.87 

5.29 

1.00 

0.80 

6.76 

12.14 

6.05 

1.51 

0.70 

4.20 

* 15.76 

7.31 

2.21 

0.60 

15.37 

9.04 

3.45 

0.50 

11.77 

15.68 

10.79 

0.40 

6.28 

* 15.73 

* 11.60 

0.30 

1.23 

8.91 

6.37 

0.20 

1.32 

The  maximum  value  of  the  Amplification  Factor 
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In  an  attempt  to  establish  the  "en"  criterion  for  the  flow  tube  system, 
t another  set  of  the  analytical  calculation  has  been  performed  for  flow 

conditions  AT  = 0°C  to  16.67°C  and  (JG  = 1.83  m/sec  to  6.71  m/sec. 

The  computed  spatial  amplification  factors  for  boundary-layer  disturbances 
are  shown  in  Figures  16  to  22.  Three  different  criteria  (namely,  ampli- 
fication ratios  of  e^,  e^,  and  e^)  have  been  used  to  estimate  where  the 
boundary  layer  would  undergo  transition.  A summary  of  these  results  is 
presented  in  Figure  23. 

It  is  noted  that  the  experimentally  quoted  transition  Reynolds  number  is 
based  on  the  assumption  that  the  boundary  layer  begins  at  a virtual  origin 
30.5  cm  upstream  of  the  test  section.  In  the  present  study,  a more  real- 
istic assumption  is  made  in  that  the  boundary  layer  is  assumed  to  begin 
at  the  upstream  side  of  the  contraction  section.  Therefore,  the  experi- 
mental data  in  Figure  23  have  been  corrected  from  those  shown  in  Figure  1 
by  this  difference  in  origin  in  order  to  compare  with  the  analytical  results. 

In  general,  the  analytical  curves  predict  the  same  trend  as  the  experimental 
data.  For  surface  overheats,  AT,  less  than  7°C,  the  e curve  appears  to 
have  a better  agreement  with  the  experimental  data  which  leads  to  the  same 
conclusion  as  the  previous  analytical  results.  Since  the  increases  in 
transition  Reynold 1 s number  in  this  flow  tube  facility  have  been  achieved  by 
increasing  the  flow  velocity  while  the  basic  geometrical  parameters  of  the 
facility  remain  the  same,  the  limitation  of  the  "maximum"  achievable  Reynold's 
number  in  the  flow  tube  may  also  be  a result  of  increasing  freestream  distur- 
bance level  as  the  flow  velocity  increases  (a  unit  Reynold's  number  effect). 
This  conjecture  needs  to  be  further  substantiated  or  reputed  in  the  future. 


It  is  also  of  interest  to  note  that  the  unstable  range  of  Tollmien-Schl ichting 
waves  with  the  flow-tube  boundary  layer  is  in  the  low  frequency  region,  ranging 
from  5Hz  to  25Hz. 


3. 05m/ sec 
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Figure  17.  Spatial  Amplification  Factors  in  Boundary  Layer  for  U =3. 05m/ sec. 
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Figure  19.  Spatial  Amplification  Factors  in  Boundary  Layer  for  U =6.41m/sec, 


Figure  22.  Spatial  Amplification  Factors  in  Boundary  Layer  for  Uo=6.71m/sec, 
AT  = 16.67°C,  T = 10°C 


Overheat 


Wall  Overheat  (°C) 

Figure  23.  Comparison  Between  Experimental  and  Analytical  Boundary-Layer  Transition  Results 
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The  dependence  of  the  transition  Reynolds  number  at  large  surface  overheat 

may  be  better  understood  by  a close  examination  of  the  spatial  amplification 

results  for  various  surface  overheat  conditions  (Figure  16  to  22).  These 

results  show  that  the  contribution  to  the  amplification  factor  in  the  unheated 

14 

section  becomes  more  important  as  AT  increases.  Based  on  the  e criterion, 
the  percentage  contribution  from  both  the  heated  and  unheated  portions  of  the 
flow  tube  are  presented  in  the  following  table. 


Surface 

Overheat 

AT  (°C) 

Percentage 
Amplification  from 

Heated  Tube 

Percentage 

Amplification  from 

Unheated  Tube 

0.00 

74 

26 

2.78 

66 

34 

5.56 

64 

36 

8.33 

34 

66 

11.11 

9 

91 

13.89 

0 

100 

16.67 

0 

100 

We  note  that  the  effectiveness  of  surface  heating  reduces  to  zero  beyond  12°C. 
To  illustrate  the  improvement  in  transition  Reynolds  number  which  might  be 
realized  if  the  extension  tube  were  heated,  calculations  were  performed  by 
assuming  the  same  surface  overheat  in  the  extension  tube.  The  computed 
amplification  ratios  for  AT  = 5.56°C  and  AT  = 8.33°C  are  shown  in  Figures  24 
and  25,  respectively.  As  shown  in  the  table  below,  at  the  end  of  the  exten- 
sion tube,  as  much  as  a 45%  decrease  in  amplification  factor  may  have  been 
achieved  when  AT  = 8.33°C.  Therefore,  a significant  improvement  of  the  ex- 
perimental performance  in  the  tube  test  may  be  expected  by  heating  the 
extension  tube. 


At  End  of  One  Extension  Tube 

X/D  = 85.5 

Surface 

Freestream 

Overheat 

Speed 

Loge(A/Ao)  with 

Loge(A/AQ)  with 

AT  (°C) 

U0(m/sec) 

Unheated  Ext. 

Heated  Ext. 

Tube 

Tube 

(with  Heated  Extension  Tube) 


Figure  24.  Spatial  Amplification  Factors  in  Boundary  Layer  for  UQ=5.80m/sec 
AT  = 5.56°C,  T = 10°C  (with  Heated  Extension  Tube) 
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4.  SUMMARY  AND  CONCLUSIONS 

This  report  documents  a complementary  analytical  study  to  the  Colorado 

flow-tube  experiment.  The  main  conclusions  are: 

I An  improved  analytical  method  for  analyzing  boundary-layer  transition 
in  the  flow  tube  has  been  developed.  The  method  includes  the  effects 
of  transverse  curvature  and  upstream  contraction  on  the  boundary-layer 
development,  coupled  with  an  iterative  technique  to  account  for  the 
displacement  thickness  correction  to  the  pressure  distribution  in  the 
flow  tube. 

I A series  of  calculations  has  been  performed  to  simulate  the  flow-tube 
experiment,  a heated  test  section  with  one  unheated  section.  At  the 
experimental  transition  location,  the  computed  amplification  factors 
vary  from  e11'6  to  e15,8  with  a mean  value  of  e14'0  which  is  much  higher 
than  the  traditional  "e^"  criterion. 

I Based  on  the  "e^"  criterion,  reasonably  favorable  comparisons  with  the 
experimental  results  are  obtained  with  this  more  accurate  computation 
scheme.  The  calculations  provide  the  explanation,  at  least  partially, 
of  the  apparent  reduced  effectiveness  of  surface  heating  at  large  AT. 

The  explanation  rests  upon  the  role  of  the  unheated  section  and  the 
effects  unique  to  a tube  boundary  layer. 

I Further  calculations  indicate  that  a significant  improvement  of  the  tube 
performance  may  be  obtained  by  heating  of  the  extension  tube.  These  re- 
sults, however,  do  not  explain  the  experimental  findings  of  Barker.  Bear 
in  mind  that  this  report  only  covers  the  analytical  study  dealing  with  main 
features  unique  to  a flow  tube.  Several  other  features  which  might  limit 
the  transition  Reynold's  number  in  a flow  tube,  such  as  the  effects  of  free 
stream  turbulence,  unsteadiness  of  the  mean  flow,  buoyancy  forces,  surface 
roughness,  surface  waviness,  suspended  particulate,  or  the  influence  of 
downstream  exit  conditions  are  not  included  in  this  study. 
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